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Abstract 

We use nonperturbative lattice techniques to study the volume-reduced "Eguchi-Kawai" version 
of four-dimensional large-A^ QCD with a single adjoint Dirac fermion. We explore the phase 
diagram of this single-site theory in the space of quark mass and gauge coupling using Wilson 
fermions for a number of colors in the range 8 < < 15. Our evidence suggests that these values 
of N are large enough to determine the nature of the phase diagram for N ^ oo. We identify the 
region in the parameter space where the {Zjy)^ center-symmetry is intact. According to previous 
theoretical work using the orbifolding paradigm, and assuming that translation invariance is not 
spontaneously broken in the infinite-volume theory, in this region volume-reduction holds: the 
single-site and infinite-volume theories become equivalent when N oo. We find strong evidence 
that this region includes both light and heavy quarks (with masses that are at the cutoff scale) , and 
our results are consistent with this region extending towards the continuum limit. We also compare 
the action density and the eigenvalue density of the overlap Dirac operator in the fundamental 
representation with those obtained in large-A pure-gauge theory. 
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I. INTRODUCTION 



The expansion of SU{N) gauge theories is a particularly useful tool for exploring 
the non-perturbative dynamics of QCD and related theories. Applications can be found 
in many areas of research, including the dynamics of confinement, issues related to the 
phase diagram of QCD, and the relation of QCD to a possible string construction. While 
large- methods offer a route to approach certain QCD-related theories with some analytic 
control, the large- iV limit of QCD itself remains unsolved. This makes lattice studies of this 
limit quite useful, and indeed these have provided a substantial body of non-perturbative 
information about large-A^ QCD. This information plays an important role in guiding and 
testing the approximate analytic approaches to large- QCD.^ 

One way to approach the large- A^ limit on the lattice is to use a straightforward gener- 
alization of the methods used to simulate QCD. The continuum and infinite volume limits 
are taken for various values of A^, and the resulting physical quantities then extrapolated 
to the large- A^ limit (typical values of A^ used are 2 < A^ < 8, but in some instances larger 
values, A^ = 10 — 16, have been used). 

In this paper we use a complementary approach commonly referred to as "large- A^ volume 
reduction". The original idea was proposed for lattice regulated theories in the seminal 
paper by Eguchi and Kawai One considers a "lattice-like" matrix model that, under 
certain assumptions, can be shown to be equivalent to a corresponding (infinite volume) 
lattice gauge theory if A^ ^ oo in both theories. By "lattice-like" we mean a model whose 
degrees of freedom take values in the group, rather than in its algebra, and that depends on 
dimensionless couplings and bare masses. The equivalent lattice gauge theory has the same 
values for these dimensionless couplings and masses. The equivalence is thus to a theory 
with a fixed cutoff, and one must take the large- A^ limit before tuning parameters to take 
the continuum limit. 

The idea of Eguchi and Kawai has spawned much interesting work demonstrating large- 
A^ equivalences between different theories. The two equivalences that motivate our present 
work are the orientifold and orbifold equivalences of Refs. . Combining these leads 

to the following result ^: 



^ References to the lattice studies can be found, for example, in Ref. 
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The large- limit of infinite- volume lattice QCD with 2Nf Dirac fermions in the 
antisymmetric representation is equivalent, for some observables, to the large- 
limit of QCD with Nf Dirac fermions in the adjoint representation defined on a 
lattice with a finite number of sites A^^. The boundary conditions in all directions 
for both fermions and gauge fields must be periodic. This equivalence holds, in 
particular, for the "single-site theory" in which A^^ = 1. 

Thus, by studying the large- A^ limit of the single-site theory with adjoint fermions, we can 
explore large- A^ QCD with fermions in the antisymmetric representation. The latter theory 
has considerable phenomenological interest because it reduces to physical QCD with 2Nf 
fermions in the fundamental representation when A^ = 3. Thus, using the equivalence above, 
we are able to study a large- A^ limit of QCD that differs from the standard 't Hooft limit. 
This alternate limit has the distinguishing feature that fermion loops are present at leading 
order in the large- A^ expansion. 

In order for the combined equivalence to hold several conditions must be fulfilled {g, 7|: 

1. The ground state of infinite- volume large- A^ QCD with Nf Dirac fermions in the 
adjoint representation must be translation invariant. 

2. The ground state of infinite volume large-A^ QCD with 2Nf Dirac fermions in the 
antisymmetric representation must be charge-conjugation invariant. 

3. The ground state of the large- A^ single-site QCD with Nf Dirac fermions in the adjoint 
representations must be (Z^)^ invariant. This symmetry is the familiar center sym- 
metry in which each of the four Polyakov loops that wind around the four Euclidean 
compactified directions is independently multiplied by a Zn factor. (Note that in the 
single-site model each Polyakov loop is composed of a single link matrix.) 

It is also necessary that all theories obey cluster decomposition; the ground state must not 
be linear combination of vacua that become disjoint at large A^. This then implies that 
multi-trace expectation values factorize at large- A^. We do not expect cluster decomposition 
to fail, but we mention this condition for completeness. 



Clearly, a crucial question is whether the three conditions listed above actually hold. 
We have no reason to suspect that the first two conditions fail, and we assume here that 
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they hold. The status of the third condition is less clear. In the following we give a brief 
summary of relevant results in the literature, from which we conclude that, while there are 
some reasons to think that the condition holds, what is needed is a direct study of this issue. 

First we note that the third condition fails when the quark masses go to infinity. In that 
limit the single-site theory becomes the Eguchi-Kawai (EK) model, for which analytical 
and numerical results show that the center symmetry is spontaneously broken at weak 
coupling 0, y, [l^. This is not necessarily a concern, however, because we are interested in 
small quark masses, and there are reasons to think that there will be a transition to a phase 
with restored center symmetry as the quark mass is lowered. In particular, in the limit of zero 
quark mass, an analysis of the continuum theory on x S^, using weak-coupling techniques 
that are valid when the radius of the 5*^ is small enough, finds that the center symmetry 
(here just Zn) is unbroken for small radius 6|. The situation at non-zero quark mass has 



been discussed in Refs. 



It appears to us that the conclusion from the last of these papers 



is that the center symmetry is broken for any non-zero mass when ^ oo. We also note 
in passing that, for very heavy quark masses, the effect of the fermions is to induce extra 
interactions between Polyakov loops wrapping around the compactified direction, and from 
this point of view, the emergent model is in the class of "deformed EK models" suggested in 
Ref. for which the center symmetry is unbroken at weak coupling for a judicious choice 
of its parameters (for a related study see Ref. 



Very recently, the calculation in Ref. p] was extended in Ref. iJ] to lattice regularization 
using Wilson fermions. The results were promising: the Zn symmetric vacuum was seen 
to have a lower energy than that of the vacua breaking Zn — > 0, for a range of lattice 
parameters that are physically relevant and that include the chiral point. The results of 
Ref. [ij] also suggest that the center symmetry may be intact even for quite heavy fermions, 
thus opening a path to study the pure-gauge theory on a lattice with one direction reduced 
to a point. We note, however, a caveat concerning the results of Ref. [l^: the possibility of 
more elaborate center-symmetry breaking to nontrivial subgroups of Zn was not considered. 
Such symmetry-breaking has been found in a different, though similar, set up, in which one 
uses a continuum regulator in the directions and a lattice regulator in the direction 15|. 
In fact, in this set up the symmetry is found to break even at zero mass. 

For completeness, we note that other approaches to large- reduction have been followed 
in the literature. Most closely related to the present work is the study of EK reduction in 
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the matrix model obtained by dimensional reduction of SU{N) supersymmetric Yang-Mills 
theory (the Nj = 1/2 case discussed below) 16|. This work differs, however, in using a non- 



16|that 
13. 



compact representation of the gauge fields. Nevertheless, the evidence found in Ref. 
reduction hold for a range of scales is encouraging. See also the related work in Ref. 

In addition, for QCD with quarks in the fundamental representation, whose dynamics 
in the N ^ oo limit are those of the pure-gauge theory, it has been found that volume 
independence does hold as long as one does not reduce the length of the box below a fixed 
physical size of 0(1 fm) [l^ (see also Ref. Other approaches to repair the center 



symmetry breaking problem of the EK model, such as the twistec 



EK models, have been found recently to fail for weak coupling 



21 



22 



20l| and quenched 



23 



24|. 



m 



As can be seen, there are no results that directly address our third condition for the single- 
site model. In principle, one could extend the perturbative calculation of Ref. [l^ to the case 
where all Euclidean directions are reduced to a point, but this would require deahng with 
infrared divergences that arise in the single-site theory. In fact, even if such a weak-coupling 
calculation were done, it would not tell us the center-symmetry realization for moderate 
couplings, where actual lattice calculations are done. For example, we might find that a 
weak-coupling calculation points to a Zjy broken phase, but that the larger fluctuations in 
the gauge fields, which occur as the coupling increases, restore the symmetry. Alternatively, 
if a one-loop calculation tells us the symmetry is intact, then a sufficiently intricate vacuum 
manifold could lead to symmetry-breaking at strong enough coupling. These possibilities 
are not just academic exercises as both were observed in related models — for example see 



the phase diagram of the model studied in Ref. [22 1. 

From the discussion above it is clear that a non-perturbative lattice Monte-Carlo analysis 
of the single-site model is required, and this is what we perform in this paper. In particular 
we focus here on the theory with Nf = 1, which is connected, through the equivalences 
mentioned above, to physical QCD with two flavors. Our results suggest that the {Zj^)^ 
symmetry is intact for a broad range of quark masses including zero. Thus we are studying 
a theory which is "within 1/A^" of the infinite volume theories appearing in conditions (1)- 
(2) above: QCD at large with one Dirac fermion in the adjoint representation, QCD 
at large A^ with two Dirac fermions in the antisymmetric representation, and, last but not 
least, QCD at A^ = 3 with two degenerate flavors in the fundamental. 
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It would also be of considerable interest to study whether reduction holds for other values 
of Nf. For Nf = 1/2, the equivalence is, in the massless limit, to the large- limit of A/" = 1 
SUSY (and to QCD with a single Dirac flavor). For two colors this SUSY theory was studied 



on the lattice by various authors 



2611 (see also the review in Ref. 



27( 1 and the work noted 



above on the reduced model in the non-compact theory [3]). The Nj = 2 case is of interest 
as a potential example of a nearly-conformal or conformal theory. Again, for a small number 
of colors, these theories have been studied on the lattice j^ . |29 |. 

The following is the outline of the paper. In Sec. [TTl we discuss QCD with Nf adjoint 
fermions — the theory that our single-site model is potentially equivalent to. We define 
the theory and describe a conjecture for its phase diagram. In Sec. Illlt we discuss the 
corresponding volume-reduced single-site theory — again defining this theory and presenting 
a conjecture for its phase diagram. Section HVl provides some technical details of our lattice 
Monte-Carlo simulations of the single-site theory. In Sec. |V] we define the observables we 
measure and use them to map the parameter space of our model, looking for regimes in 
which the center symmetry is intact. Section |Vl] includes a restricted set of results of 
physical interest, such as certain eigenvalue densities of Dirac operators. We summarize our 
study in Section IVIII and discuss possible future directions of study. 



II. LATTICE QCD WITH ADJOINT FERMIONS 

If reduction holds, the single-site matrix model discussed in the next section is equivalent, 
at large A^, to lattice QCD with fermions in the adjoint representation in (arbitrarily) large 
volumes. Here we discuss the properties of the latter theory, so that we know what we 
should expect to find in the single-site model if reduction holds. 

The orbifold construction that underlies this potential equivalence implies that the form 
of the lattice actions is the same in both the reduced and unreduced theories. We use the 
Wilson gauge action and Wilson fermions for the matrix model, and so discuss below the 
same action for the large-volume theory. We thus consider a gauge theory in four Euclidean 
dimensions, with lattice spacing a and sites, and whose path integral 

^adj = / DUDi^Di, exp (^gauge + ^D^i,). (2.1) 
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The gauge action is 

-5gauge = 27V6 5^ReIVC/p, (2.2) 

p 

where P labels plaquettes, Up is product of SU {N) link variables around the plaquette, and 
b is the inverse 't Hooft couphng that is kept fixed as is increased 



{g'N)-' . (2.3) 



The Nf Dirac fields ip carry implicit spatial, spinor and adjoint color indices. We use the 
lattice Wilson- Dirac operator 



X] (1 - 7m) U^,^i ^y,x+n + (1 + 7m) Sy 



(2.4) 



where x and y label sites, and k is the usual hopping parameter, related to the bare quark 
mass by 

K = . (2.5) 

8 + 2amo ^ ' 

The boundary conditions on both gauge and fermion fields are taken to be periodic. 

In a literal implementation, C/^^ would be the adjoint representatives of the SU{N) 
matrices U^^fj, appearing in the gauge action. Thus would be a matrix of dimension 
(A^^ — 1), and i/j an (A?"^ — l)-dimensional color vector. We find it simpler to place the 
fermions in the reducible N i^i N — adj. © 1 representation, so that they have A^^ color 
components, and are acted on by A^^ x A^^ matrices. As will be seen shortly, the additional 
singlet component decouples from the dynamics. 

We denote fundamental representation color indices by lower-case letters, e.g. a, 6 e 
[1, A/"], and N i^i N indices by upper-case letters, e.g. A,B e [1, N'^]. We choose a basis in 
which the latter index is composite: 

A={ai;a2), (2.6) 

and, correspondingly, in which 

This acts on a fermion field with indices ipB = ^(fei;62)- By change of basis we can bring U'^ 
into block diagonal form, with a 1-d block for the singlet, and an (A^^ — l)-d block for the 
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desired adjoint part. The singlet part is, however, unity for all U E SU{N). This implies 
that 

det DwiU^^j^) = det L)vi/(f/adj.) x constant , (2.8) 

showing that we obtain the path integral of the desired adjoint theory up to an irrelevant 
overall constant. To calculate fermionic expectation values one would, however, need to 
remove the singlet component. 

At tree-level in perturbation theory, the quarks become massless when /t = 1/8, a result 
which holds for any representation. At higher order, and non-perturbatively, amo is addi- 
tively renormalized, and the true chiral point occurs at Kc > 1/8- This is because Wilson 
fermions break chiral symmetry explicitly for any nonzero lattice spacing. The physical mass 
is then given by 

= V - 2^) ' ^2.9) 
with a multiplicative renormalization factor that depends on the scheme chosen to define 
the mass. This factor is finite and of 0(1) as long as we consider lattice spacings that are 
not too small, as we do in practice. 

QCD with a single adjoint Dirac fermion is asymptotically free (and quite far from the 



conformal window explored in Refs. 28|, |29|), and is expected to confine and spontaneously 



break its chiral symmetry. Let us first consider the theory with a chirally symmetric regula- 
tor. Because the fermion representation is real, its global symme try gro up in the chiral limit 
is SU{2Nf) and not SU{Nf)L x SU{Nf) (for example, see Refs. [30|, l31|)- This symmetry is 
conjectured to break spontaneously to the flavor group S0{2Nf), generating 2A^| + Nj — 1 
Nambu-Goldstone bosons. For Nj = 1 we therefore expect two massless modes. The lattice 
theory, with Wilson fermions, respects only the S0{2Nf) flavor symmetry (which can be 
most easily seen by rewriting the action using Majorana fermions). Chiral symmetry is 
restored only in the continuum limit, and requires tuning k appropriately. We discuss the 
situation at finite lattice spacing below when we sketch the phase diagram. 

The theory also has a {Z^Y center symmetry, each factor corresponding to multiplying 
all the link elements U^^^o on a slice of fixed coordinate Xq in the /iq direction by an element 
of Zn- Since the fermions are in the adjoint representation, they are neutral under this 
symmetry. We expect this symmetry to be unbroken for large volumes, as it is in the 
pure-gauge theory. The key issue is whether it remains unbroken in small volumes. As 



7 



discussed in the Introduction, analytical arguments in the case of a single short direction 
indicate that it will remain unbroken at small volumes if the boundary conditions on all 
fields are periodic. This is in contrast to what happens when the fermions have antiperiodic 
boundary conditions, when we do expect the center symmetry to be broken when one passes 
through the finite-temperature deconfinement transition as the physical length of one of the 
directions is reduced. 



A. Phase diagram of adjoint lattice QCD with Nj = 1. 



As far as we know, there have been no lattice studies of gauge theories with Nf = 1 adjoint 
fermions, even for = 2 or 3. There has been extensive work for Nf = 2, with N = 2 and 



3, both for non-zero temperature (see, e.g., Refs. [32|, |33|) and at zero temperature 28l. |29|. 
As discussed in the Introduction, the latter theories are expected to be conformal or nearly 
conformal, and are not likely to provide useful guidance for the Nf = 1 case. There has 
also been much work on the Nf = 1/2, N = 2 theory with both Wilson fermions (as 



reviewed in Ref. 
fermions 



25 



3J]; see Ref. 35|] for recent progress) and, more recently. Domain- Wall 



261 . Here, again, the results do not obviously apply to Nf = 1 theories, 



because the target Nf = 1/2 theory is, in the chiral limit, supersymmetric. 

In fact, it may well be that QCD — physical QCD with fermions in the fundamental 
representation — is the theory whose dynamics is most similar to that of the Nf = 1 adjoint 
theory. Assuming so, we make the educated guess for the phase diagram shown in Fig. [H 
Although we are particularly interested in the large- A^ limit, we expect this sketch to hold 
also for small values of A^. 

Let us explain the features of this diagram. The solid (blue) line labeled Hc{b), and ending 
at (k, b) = (0.125, oo) is the "critical" line (or possibly lines, as discussed below) along which 
to the fermions attain their minimal mass. The continuum limit in which there are light 
fermions is obtained by approaching the end point of this line (in a way which depends upon 
the desired fermion masses). 

Close enough to the continuum limit (i.e. for large enough b), the theory in the vicin- 
ity of the critical line can be analyzed using chiral perturbation theory including lattice 



artifacts 
SU{2) - 



36(1 . The analysis is similar to that for QCD, but must be generalized to the 
50(2) chiral-symmetry breaking pattern. In fact, the corresponding case for 
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0.125 



OO 








-0.36 - 



0.25 



FIG. 1: Conjectured phase diagram in the k — b plane for a QCD-hke theory with a single Dirac 
adjoint fermion. Details are discussed in the text. 

Nf = 2 adjoints [SU{A) S'0(4)] has been worked out in Ref. [sT], and the generahzation 
to our case is straightforward. One finds that, as in QCD, there are two possibihties: ei- 
ther there is a hne of first-order transitions, at which the two degenerate pseudo-Goldstone 
"pions" attain their minimal (non-zero) mass, or there are two lines of second-order tran- 
sitions, at which the masses of both pions vanish. Between these second-order lines lies an 
Aoki-phase 38|, in which the 5*0(2) flavor symmetry of the lattice theory is spontaneously 
broken, so that one of the pions is exactly massless, while the other is massive. As 6 ^ oo, 
the non-zero pion masses in both scenarios go to zero linearly with a. In addition, the width 
of the Aoki-phase, if present, shrinks as a^. 

One cannot predict which scenario applies without a non-perturbative calculation, and 
the result depends on the details of the action. The numerical results discussed below suggest 
that the transition is first order (and thus that Kc is a single line) for b > 0.1 (which is the 
smallest value we use). It is worth stressing that, in either scenario, one can approach the 
continuum limit on either side of the transition(s), i.e. with k < Kc(min) or k > Kc(max). 
The two sides have identical long-distance physics as a — 0. 

As one moves to stronger coupling, terms of higher order in a become important, and it 
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is possible that one changes from a first-order transition to having an Aoki phase, or more 
exotic possibihties. Once at very strong couphng, g'^N ^ 1, one can show that there will 
be an Aoki phase, and that 1/4 when 6 — > 0. This is because the analysis of lattice 

QCD in this limit carried out in Ref. holds also for the theory with fermions in the 
adjoint representation.^ Because of this, we have shown a region (solid [blue] shading) of 
Aoki-phase for 6 < 0.1. 

For > \ there will be additional phase structure (the Aoki phase "fingers"), with 
additional critical lines along which the continuum theory has more than one light Dirac 
fermions (four or six, depending along which critical line one takes the continuum). These 
continuum theories are not asymptotically free, and are not interesting for our purposes of 
connecting to physical QCD. Thus we have restricted our attention to k < i. 

The other feature shown in Fig. [1] is the approximately horizontal (red) dashed line 
punctuated with question marks. This indicates a possible "bulk" transition, so we label it 
&buik('t)- For > 5 such a transition is known to be present at k = (i.e. the pure-gauge 
theory) and to be strongly first-order (for lower values of A^ it is a crossover). We therefore 
expect that it persists as a first-order transition line for some distance out into the k, — b 
plane. We emphasize that this transition is a lattice artifact and so the bare coupling at 
which it takes place approaches a non-zero (and non-infinite) k— dependent value for infinite 
volume and infinite A^. For example, it occurs for k = at 6buik('^ = 0) ^ 0.36 as indicated 
in the Figure.^ This transition is also not associated with any symmetry breaking and so can 
end anywhere in the (k — &) plane. In the absence of any data we do not know what happens 
and so decorate the dashed line by question marks. In any event, since the continuum limit is 
at 6 = oo, lattice simulations aiming to approach that limit should be made for b > 6buik('t)- 



^ This holds despite the fact that the symmetry breaking pattern is different. For example, one of the pions 
which becomes massless is created by the "diquark" operator ip'^Cj^tp. For g'^N ^ 1, the propagator for 
this diquark has, in a hopping-parameter expansion, both quarks hopping along the same path, so that 
gauge matrices completely cancel from the "pion" propagator, just as for quark-antiquark pair in QCD in 
this limit. One can furthermore show that the Dirac-matrix factors are the same in the two calculations. 
Thus the propagators in the two theories are proportional, and so the corresponding "pion" masses vanish 
at the same k. 
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The value of t'buik('^ = 0) was measured, for example, by hysteresis scans [40|. A result is also given in 
Ref. 4l|, but using the simulations of the "twisted EK" model, a variation of the EK model that was 



for the values TV < 64 used in Ref. 



recently invalidated for large enough values of N [23,l23|. Nonetheless, it is possible that reduction holds 



4l|, and that the estimate for 6buik('v — 0) obtained there is reliable. 
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III. THE VOLUME REDUCED THEORY AND ITS PHASE DIAGRAM 



We now turn to the single-site theory which is the focus of the present work. It is defined 
simply as the L = 1 version of the construction in Eq. fl2.1l) . and is the generalization of 
the original EK model obtained by adding adjoint fermions. The degrees of freedom in the 
model are the four SU{N) matrices t/^=i_4, and the 4A^^-component Grassmann variables 
ip and ip. The action of this adjoint EK (AEK) model is 

^AEK = S^K + ^D'w^, (3.1) 

Saek = 2NbJ2 ReTr Ul Ul , (3.2) 

^l<u 

D'^^ = l-'^E[(l-7.)t/i'+(l + 7.)f/;''] • (3.3) 

For the purpose of Monte-Carlo simulations, we formally integrate over the fermions and 
evaluate expectation values with the following path integral 

Zaek = I i{DU, exp [^EK + log det (D]^)] . (3.4) 
JSU{N) ^J-^ L V yj 

The relevant symmetries of Zaek are the same as those of the original EK model. They 
are the remnant of the gauge symmetry 

V/i : U^^VtU^n^ with VteSU{N), (3.5) 

as well as center transformations applied independently to the four link matrices 

U^z""^ with z = e^'^'l^ and G . (3.6) 

As explained in Sec. [H large-A^ equivalence holds as long as the iZ^y" symmetry in 
Eq. (13. 6p is unbroken. We recall here how this equivalence works in detail. This equivalence 
states that appropriate expectation values in the reduced theory become identical when N — >■ 
oo to those in the infinite- volume theory defined by Eq. (12.11) . Appropriate expectation values 
in the large volume theory are the connected correlators of (Z/vr)^-invariant and translation- 
invariant operators. These are mapped into operators in the reduced theory following the 

6|. For example, consider the large- volume expectation value of the 



prescription of Refs. 0, 6 
plaquette. 



«-^]^E(t''fp)-..- (3-7) 
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The notation (, )^^^. means that we calculate expectation values in the ensemble defined by 
the partition function in Eq. (12. ip . Np is the number of plaquettes, which in four dimensions 
is equal to GL"^. The corresponding single-site expectation value is 

^^lEi^'U.U^Ul Ul )z,,, (3.8) 

so that, in fact, 

Ured. = u{L = 1) . (3.9) 
The meaning of volume reduction is that 

m(6, k) = Mred.(^5 '^) (3.10) 

when — > oo in both theories. 

Our aim in this paper is to find the regions of the 6 — k plane in which the ground-state 
of the single-site model is invariant under the (^at)^ center symmetry, so that equivalences 
of the form of Eq. (13.101) hold. We first collect what is known about the single-site theory, 
together with some conjectures, into a phase diagram. 

For infinitely massive fermions (i.e. for n = 0), our theory becomes the original EK 



model. This is known to break the (Z^)^ symmetry 
evidence suggests that the transition is first order [sl, [iC 



'or b > bv.K ~ 0.19 and numerical 



421]. A crucial issue is then whether. 



for b > 6ek; an increase in k can lead to the restoration of the center symmetry. This is 



]J] studied 



what one would expect based on the results of Refs. 0, Q . Specifically, Ref . 
a lattice theory similar to Eq. (13. ip but in which only one of the Euclidean dimensions is 
reduced to a point. This weak-coupling calculation found that, for the Nf = 1 theory, the 
center symmetry is broken for k = 0, but restored once k, ^ 0.04. There is, in addition, a 
small intermediate phase between the two in which the symmetry is broken down to a 
Z2 subgroup. Finally, when k grows even more, to values above ~ 1.4, the Z^v symmetry is 
completely broken again. 

It is unclear how the fact that in the current paper we study a theory where all lattice 



directions are reduced to a point changes the results of Ref. ij]. Nevertheless, assuming 



that the results of Ref. [14] provide a qualitative guide, we are led to conjecture the phase 



diagram shown in Fig. O In the following we describe the features of this diagram. 

First, consider the solid (black) line which begins at 6ek ~ 0.19 on the 6-axis, bends up, 
and ends on the top of the diagram. Our conjecture is that this separates a center-symmetry 
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0.125 



oo 



-0.19 








0.25 



FIG. 2: Conjectured phase diagram for the single-site model in the large-A^ limit. Details are 
discussed in the text. 



broken phase to the left (shown shaded) from an unbroken phase to the right. This is based 
on what we know about the EK model, and the assumption that the results of Ref. [l^ 
apply qualitatively to our model. (Note that there can only be a phase transition when 
oo, while for finite N this will be smoothed out into a crossover.) It is possible that 
this line actually has a non-zero width and contains intermediate partially broken phases.^ 
We stress that we do not know the value of k at which this line hits the b = oo axis, but our 
conjecture is that this value is smaller than (likely significantly smaller than) k = 0.125. 

Towards the right-hand side of the plot there is another solid (black) line which also 
separates the conjectured central center-symmetric phase from a shaded symmetry-broken 
phase or phases. Here we actually expect a very complicated phase diagram, based on the 
fact that there are regions where the corresponding infinite volume theory has multiple light 
fermions (as discussed above), and the results of Ref. [3]. Since this is not our region of 
interest, we have not tried to fill in the details. What is important here is the conjecture that 
this symmetry-broken region lies (significantly) to the right of the critical line (consistent 



For example, this was seen in Ref. ^] at k — 0. 
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with the results for Ref. 




Given these conjectures, we are left with a central "funnel" in which reduction holds, and 
thus have reproduced both the critical "line" Kc{b) and the possible bulk transition line from 
our conjecture for the infinite-volume theory. Fig. [T] The key question for our numerical 
investigation is whether this central funnel, and in particular its upper part (where one can 
take the continuum limit) is actually present. We stress again that the precise position of 
the boundaries of this funnel are supposed to represent only the conjecture that there is a 
generous region on either side of Kc in the symmetry-unbroken phase. 

Reduction does not hold in the regions where the (Zn)'^ is broken. For this reason we 
have changed the character of the dashed (red) &buik(^) line outside of the central funnel. In 
particular, we are not aware of numerical evidence for a bulk transition for k = beyond 
6ek and did not ourselves observe one. Thus, we end the line away from the 6-axis. We 
stress that it is not a priori known whether there is a bulk transition at all for any value of 

K. 

In the next sections we study the theory defined by Eq. (13. 4p using non-perturbative 
Monte-Carlo simulations, and indeed find a phase diagram similar to that appearing in 
Fig.H 



IV. NON-PERTURBATIVE LATTICE STUDIES: TECHNICAL DETAILS 



We study the path integral in Eq. (13.40 using Monte Carlo simulations. The weight 
function is 

P(U) = e^^K(^) det D'^'^iU) , (4.1) 

which is integrated using the SU{N) Haar measure for each link: n^=i dU^. For a single 
Dirac fermion in the adjoint representation, det D^^ is real and positive, so that P{U) 
can be treated as a probability density. The reality of the determinant follows as usual 
from 75 hermiticity (75-Dvi/*^75 = (D^'Y). Positivity follows because the fermion is in 
a real representation, which allows the action to be rewritten in terms of two Majorana 
fermions, each of which gives a Pfaffian when integrated out. The Pfaffian is real, though of 



indeterminate sign, but its square is necessarily positive 



34|. 



To produce the field configurations we use a standard Metropolis algorithm. Following 
[ I 

Cabibbo and Marinari ^], our proposed changes are obtained by multiplying the links by 
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matrices living in SU{2) subgroups of SU{N). For each subgroup, we propose five changes, 
and then run through the N{N — l)/2 SU{2) subgroups in turn. We repeat this for each of 
the four hnks. To calculate the change in P{U) we simply calculate the determinant anew 
after each suggested change in the links. The 10A^(A^ — 1) proposed changes just described 
constitute what we call a "model update". We perform measurements every five model 
updates, after a number of initial "thermalization" updates. We found that we could attain 
acceptance rates of 50-60%.^ 

Since the cost of calculating the determinant of an M x M matrix scales like M^, the 
cost of each of the SU{2) updates scales like A^^, and the overall cost of a model update 
scales like A^^. This means that, for a fixed number of model updates, a calculation with 

= 15 is 25 times more expensive than one with N = 10. Our resources for this calculation 
were very modest — roughly an average of four Intel(R) 6700 @ 2.66GHz CPUs. We did not 
attempt to parallelize our code, and since we have a single lattice site, it is not clear to 
us if this is possible. On a single CPU, an SU (10) calculation including 50 thermalization 
model updates and 100 measurements (550 model updates altogether) took around 3| hours. 
We also found that the asymptotic scaling held to reasonable approximation (so that 
gathering 500 model updates in SU{15) takes about three days, explaining the modest data- 
set we collected for that gauge group). We note that it may not be necessary to update 
all the SU{2) subgroups of SU{N) in a single model update — such a procedure is perhaps 
more suitable for pure-gauge models whose computational scaling law is a moderate N^. 
Nevertheless, because our calculation is first of its kind, we aimed to be conservative and to 
avoid auto-correlations between successive measurements to the extent possible. 

In the following, all errors have been calculated using the jackknife procedure. In some 
measurements, we varied the bin size and chose a value for which the statistical error satu- 
rated. In other measurements, however, we worked with a fixed bin size, and in these cases, 
based on our experience with variable bin sizes, we multiplied the resulting statistical error 
by a factor of two. This factor is chosen to be conservative. 

^ We did so by preparing a list of 800 random SU{2) matrices (and their inverses) such that their traces 
were Gaussianly distributed around the unit matrix with a width w that decreased with increasing b and 
N. For example, for N = 8, w was ^ 0.5 for h = 0.1, ^ 0.05 for b = 0.30 and ~ 0.01 for b = 1.0. Compared 
to SU{8), the width w of SU{15), was decreases by about 20%. The distribution in the other, angular, 
directions in the SU (2) manifold were uniform. 
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Gauge group 


b 


K 




0.35 


0.065 - 0.22 


SU{8) 


0.10 - 0.50 
1.00 


0.001 - 0.40 
0.05 - 0.20 




10"^ - 1.00 


0.03,0.04 




0.30 


0.01-0.20 




0.35 


0.065 - 0.22 


SU{10) 


0.40 


0.13 - 0.18 




0.40 


0.13-0.18 




0.50 


10-^ - 0.495 




0.30 


0.03,0.06 


SU(11) 


0.50 


0.03 - 0.33 




1.00 


0.09 


SU{13) 


0.50 
0.35 


0.06 - 0.155 
0.15-0.20 


SU{15) 


0.50 
1.00 


0.06 - 0.155 
0.09 



TABLE I: Details of runs that consisted of 100-400 measurements. 

Gauge group b k Number of configurations 

0.35 0.1275,0.150,0.155 
SU{10) 0.50 0.1275 1000 

1.00 0.09,0.1275 
SU{13) 1.00 0.09 3700 

TABLE II: Details of longer runs. 

We conclude this section by presenting in Tables [T] and [Tll the details of the gauge con- 
figurations we have accumulated. We typically used 50 model updates for thermalization 
before starting measurements. 
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V. NON-PERTURBATIVE LATTICE STUDIES: RESULTS 



In this section we present the results of our Monte-Carlo studies and focus on the way the 
center symmetry is realized in various regions of phase diagram. The results in Sees. IV BtIV El 
are from the runs in Table [H while those in Sec. IV Fl are from the longer runs in Table [ITl 

A. Definition of observables 

The observables we use to map the phase diagram are listed below. They were chosen to 
probe the large number of possible breaking patterns of the center symmetry (as described, 



for example, in Ref. [l^), as well as to detect phase transitions that do not involve a change 
in the realization of the center symmetry. 

1. Plaquette defined in Eq. (13.81) . 

This observable is not an order-parameter for center symmetry, but allows us to detect 
transitions that are unrelated to the center realization. These include a possible bulk 
transition at 6buik('^); which would separate the lattice strongly-coupled phase from 
the continuum one (see Fig. [2]), and the Kc{b) line to which one needs to tune to obtain 
the minimal quark mass. 

2. The expectation values of the traces of four link variables = -^trt/^, and their 
magnitudes |P^|. We often refer to these variables as "Polyakov loops". 

The are order parameters for the complete breaking of the center symmetry, though 
they are not sensitive to partial breaking. The |P^| are probes of large- factorization: 
they scale like 1/iV if factorization holds, but like if it breaks down. 

3. The expectation values of the traces of the twelve 'corner variables' M^^, = j^ti UfJJy 
and M^_i, = -^tiU^Ul (with /i > z/), and their magnitudes |M^^-|-i,|. 

These observables were identified in Ref. as probing a nontrivial form of symmetry- 
breaking characterized by correlations between the link matrices, U^, in different di- 
rections. 

4. For some parameters we also calculate the following set of {L + 1)^ traces 

ir^ = Itr f/f^[/f f/f, withn^ = 0,±1,±2,...,±L, (5.1) 
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FIG. 3: Scatter plots of the four Polyakov loops for = 10 and k = O.The values of the coupling 
are b = 10"^ (left), b = 0.3421 (middle) and b = 0.5 (right). 

where f/^" = f/^". We take L = 5 and so calculate 14, 641 different averages for each 
gauge configuration. 

Like the M^^±^, these traces are order parameters for intricate symmetry breakdown 
schemes [such as {Z^Y or {ZjqY (^tv)^ ^ etc]. 




B. Results for k = 0: infinitely heavy quarks 

We begin by making a connection with the original EK model (obtained when k = 0), 
for which the behavior is known from previous work. This corresponds to moving up the 
left-hand axis of Fig. O We show in Fig. [3] scatter plots of the four Polyakov loops for three 
values of 6 at = 10. The smallest value, h = 10~^, is clearly in the strong-coupling regime 
where, as discussed in Sec. IIIIl we expect the center symmetry to be intact. That this is 
the case is shown by the clustering of around the origin. As b increases, the distribution 
spreads out while remaining centered on the origin (not shown), until one reaches b ^ 0.19, 
at which point center-symmetry breaking is observed.^ An illustration of the behavior well 
inside the symmetry-broken phase is shown in the right two panels, which are for b = 0.3421 
and 0.5. The Polyakov loops are seen to mainly fluctuate around elements of the center of 
SU{10), up to an overall scaling factor, i.e. {Pti)zAEK ~ e2«n^7r/io^ with po ^ 0.7. For 
b = 0.3421 there are also tunneling transitions between different center phases. Although 
at finite one does not have a phase transition in the single-site model, these figures show 



As noted above, this occurs by a cascade of transitions, analogous to those observed in Refs. 42l |. 
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that one can nevertheless observe the putative phase structure even at moderate values of 
A^. Indeed, we obtained similar results (not shown) for = 8. 

The nature of the transition is shown in more detail in Fig. HJ which displays the results of 
a scan in 6 at k = for = 10/ There is a rapid rise in the plaquette starting at b ^ 0.19, 
and a corresponding increase in the average magnitude of the Polyakov loop. Note that this 
"transition" is not related to the "bulk" phase transition of lattice gauge theories, for the 
latter does not break the center symmetry. We find that, while = 10 is large enough to 
show a clear indication of the phase transition, it is too small to see a true hysteresis curve 
(i.e. with a nonzero range of metastability). 
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FIG. 4: Plot of the average plaquette ([blue] pluses) and the magnitude of the Polyakov loop P^=i 
([red] crosses) at k = as a function of 6 for = 10. 



^ In this and subsequent scans, the number of thermalization runs at each value of b was 50, and a 100 
measurements were performed. The initial gauge configuration for each value of b was the final gauge 
configuration at the preceding value of b. The first b simulated was that at 6 = 0, where the initial 
configuration had ?7^=i_4 = 1. 
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C. Results for k > 0: dynamical quarks 



We now ask what is the effect of increasing k from zero, i.e. of decreasing the mass of 
the adjoint quarks from infinity. Most interesting is to measure this effect for values of b 
for which the center symmetry is broken at k = 0. As an illustration we show in Fig. O 
the scatter of in the complex plane for a range of k values at 6 = 0.3 and = 10. We 
observe that symmetry remains broken at k = 0.03, but appears to be restored for k = 0.06 
and above. That the symmetry is restored at such a small value of k is consistent with the 
weak-coupling analysis of Ref . [l^ . It is encouraging for our purposes since it indicates that 
even very heavy adjoint quarks are sufficient to recover the reduction that is absent in the 
EK model. 

We find that the restoration of the center-symmetry also holds both for weaker 
couplings — as in the b = 0.5 data shown in Fig. [6]-and for larger values of N — as exem- 
plified by the = 15, 6 = 0.5 results shown in Fig. [71 (Note that the first non-zero value of 
K in the A^ = 15 plots is 0.06 and not 0.03, unlike in the previous plots.) We observe that 
for K, > 0.06 the fluctuations in the Polyakov loops do decrease as A^ increases, qualitatively 
consistent with the expected 1/A^ fall-off if one has large- A^ factorization. Another way to 
see this is to look at the Monte-Carlo time history of the Polyakov loops. We present ex- 
amples for A^ = 8, 10 and 15 at 6 = 0.5 in Fig. [HI We observe that the fluctuate around 
zero, with an amplitude that decreases as A^ increases.^ 

So far, our results are consistent with the left-hand part of the conjectured phase diagram 
of Fig. [21 i.e. with the center-symmetry broken phase ending for small k. To investigate 
further, we have done scans in b for fixed non-zero n. The left panel of Fig. M compares, for 
A^ = 10, the scan of the plaquette for n = (already shown above) to that for k, = 0.0925 
(well into the putative symmetry- restored phase). The right panel shows the corresponding 
scans for |P^=i|. We see that the would-be center-breaking transition at 6 ~ 0.19 for k = 
is replaced by a new structure at larger values of 6 ^ 0.28 — 0.30, and that this new structure 
is not accompanied by an increase in |Pi|. Thus, although there may be a transition for 
K = 0.0925, it is not associated with center-symmetry breaking. It is in fact not clear from 
our study whether there is such a "bulk" transition at all. In either case, it is clearly safer 

® The relatively long decorrelation times evident in these plots suggest that longer runs may be needed for 
some parameters. We return to this point in Sec. IV Fl 
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FIG. 5: Scatter plots of the four Polyakov loops for b = 0.3 and N = 10, The value of k increases 
from the top-left plot to the bottom-right plot in the order 0.01, 0.03, 0.06, 0.11, 0.165, and 0.185. 
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FIG. 6: As in Fig. [S] but for b = 0.5, and for k values 0.0001, 0.03, 0.06, 0.12, 0.1475, and 0.16 
(from top- left to bottom-right). 
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FIG. 7: As in Fig.[6]but for N = 15 at b = 0.5. The values of k are 0, 0.06, 0.09, 0.1275, 0.145 and 
0.155 (from top-left to bottom-right). 



to work on the weak-coupling side of this possible transition when trying to make contact 
with the continuum. This does not, however, appear to pose any difficulty, since our data 
is consistent with the Zj^j symmetry being intact even at 6 = 0.5 (and, less convincingly, at 
6 = 1.0 as well — see below). 



D. Looking for the chiral point 

We now turn our attention to finding the critical line (or lines), Kc{b), that were discussed 
in Sees. [Ill and IIIII and appear in both Figs [1] and [2l To do so we performed scans for 
^ K < 0.5 at several values of b and A^. We begin by discussing the results of a hysteresis 
scan (i.e. a scan done by both increasing and decreasing in n) for SU{8) at 6 = 0.35. The 
results are presented in Fig. [TOl The left panel of the figure shows the plaquette and reveals 
a very strong hysteresis, indicative of a first order transition. The average magnitude of the 
Polyakov loops, shown in the right panel, displays a much weaker hysteresis, with (|-P^|) ^ 1 
on both sides of the transition. Based on this, and on the behavior of the scatter plots of 
the and M^i, (not shown), we conclude that this transition does not involve any breaking 
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FIG. 8: Monte-Carlo time history of the Polyakov loop P^=i for = 8, 10 (upper panels) and 
= 15 (lower panel) at 6 = 0.5. The real part of Pi is shown by [red] pluses and the imaginary 
by [blue] bursts. The values of k are 0.148 in 5C/(8), and 0.155 in SU{W) and SU{15). 

of the center-symmetry. 

We next show, in Fig. [TT|, a compilation of all our = 8 results for b = 0.1 — 1.0. We 
observe that, as b increases, the transition shifts to smaller values of k, and the discontinuity 
in the plaquette decreases. These two features are qualitatively consistent with expectations 
for the critical line K,c{b) in the first-order scenario discussed in Sees. [ITl and IIIII In particular, 
we expect that the transition should interpolate between k, = 0.125 at 6 = oo and k = 0.25 
at 6 = 0, and our results are consistent with this expectation.^ We also expect that the 

^ If there is an Aoki-phase for small b this transition line would correspond to the position of the left-hand 
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FIG. 9: Scans of the plaquette (left panel) and the magnitude of the Polyakov loop P^=i (right 
panel) for = 10 at k = ([red] pluses) and k = 0.0925 ([blue] crosses). 



transition should weaken rapidly as b increases, as observed. 

It is also important to investigate the dependence. If, when N ^ oo, there is a first- 
order transition, we would expect an increasingly wide region of metastability as increases. 
We have studied this at 6 = 0.5 for N = 8, 10, 11, 13 and 15, with results shown in Fig. [121 
We see only a very weak dependence on away from the transition (indicating that our 
results for the plaquette are close to their N = oo values), while there is some evidence for 
increasing metastability as A^ increases. This is most notable for A^ = 10 ([black] crosses) 
compared to A^ = 8 ([red] pluses). 

If our interpretation is correct, then the "pions" composed of adjoint quarks should have 
a minimal, non-zero mass along the critical line, and the long distance physics on both 
sides of the transition should be the same. It is beyond the scope (and resources) of the 
present calculation to test these claims directly. Clearly this is an important issue for further 
work. In this regard, it is useful to convert our values of b into the corresponding values of 
/3 in a standard SU (3) simulation with Wilson gauge and fermion action. The relation is 
/3 = 2N'^b, so that b = 0.35, 0.5 and 1 convert to P = 6.3, 9 and 18, respectively. In pure 



boundary of this phase. We have not done detailed studies at small b to elucidate a possible region of 
Aoki-phase. 
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FIG. 10: Hysteresis scans in k for = 8 at 6 = 0.35, for the average plaquette (left panel) and 
the average magnitude of the Polyakov loops (right panel). Scans with increasing(decreasing) k 
are shown using [red] pluses ([blue] crosses). 
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gauge simulations, one enters the weak-coupling regime at /5 ~ 6, and this crossover value 
is reduced by the presence of fermions. Consequently, our calculations at 6 > 0.4 are well 
inside the weak coupling region. 

It is interesting to compare the phase diagram we have found thus far to that observed 



for the 



are 



ge- volume Nf = 2, N = 2 theory in Ref. [29j]. In scans of the plaquette versus 



K, Ref. 291] shows a similar first-order transition for stronger couplings, but also finds that 
the line of such transitions ending at the point 6 = /5/8 ~ 0.25, k, ^ 0.19. This is a very 
different behavior from that we have observed, consistent with the underlying physics being 
itself quite different. 



E. Exploring larger k values 

We have also performed some scans at larger values of n. Here we are outside the regime 
which is connected by reduction to physical QCD (which is roughly 0.05 < k < 0.2), but this 
region is interesting for several reasons. First, we want to determine the boundaries of the 
region in which the center-symmetry is unbroken, so that reduction holds. We also want 
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to make a connection with the one- loop computation of [14!], which found a Zjy-symmetry 
breaking transition for k > 1.4. Finally, it is simply interesting in its own right to understand 
the phase structure of this single-site model. 

What we find is a complicated set of phase transitions, involving partial breaking of the 
center-symmetry. These are reminiscent of the transitions found in the one-loop potential 



studies of Ref. [iSj. We have not fully untangled the phase structure, and will thus only 
present a sampling of our results along with some conjectured interpretations. 

As an example, in Fig. [13] we present scatter plots of the observables (defined in 

Sec. IV Ap . for = 10 and b = 0.5. The value of k varies from 0.0001 (where the theory is 
essentially pure-gauge), through the intermediate values of k discussed in Sec. IV Dt up to 
large values (by which we mean n > 0.2). The corresponding Polyakov loop scatter plots 
are presented in Fig. [T3] (note that there is some overlap with Fig. [6]). We observe that, 
for K, < 0.2, the M^^±y behave similarly to the Polyakov loops: they show center-symmetry 
breaking for k ~ [if the links, U^, have phases close to e^*"'''^/^ then the M^^±^ have 
phases close to e^*("f ^"'')'^/^], but then fluctuate around zero for intermediate values of k, as 
exemplified by the results shown for k, = 0.1275. The fluctuations are, however, much larger 
for the M^i, than for the Polyakov loops. 

A new behavior is seen for n > 0.2. As illustrated by the results at k = 0.245 and 0.275, 
the Mf^i, show symmetry breaking (although with tunnelings between "vacua") while the 

do not. This indicates a mode of center-symmetry breaking involving the correlation of 
eigenvalues of links in different directions, while the eigenvalues themselves remain uniformly 
distributed. This is the pattern we observed in the quenched EK model {2^, and shows the 
importance of using order parameters other than the Polyakov loops. 

The onset of this behavior is examined in more detail [now for SU{11)] in Figs. [15] and 
[T6l The flrst two panels show the distribution of M^^, spreading out, with some of the them 
having an almost flxed phase, while the remain close to the origin and show no signs 
of symmetry breaking. The bottom two panels in each flgure show that, as k is increased 
further, all the M^^^ show a clear symmetry-breaking pattern, while the P^ start to move 
away from the origin. This is also seen in the penultimate panels of Figs. [13] and [TU 

By studying the distribution of the individual P^ and M^^, we have determined a possible 
explanation for some of these scatter plots in terms of the behavior of the eigenvalues of the 
link matrices. As an example, consider the case of = 10, 6 = 0.5 and /t = 0.275, shown 
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FIG. 13: Scatter plots of the twelve M^,±^ for iV = 10, 6 = 0.5. Here k = 0.0001, 0.1275, 0.245, 
0.275, 0.29, and 0.495 (running from top-left to bottom-right). 
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FIG. 14: Scatter plots of the four for the same data set as in Fig. [131 
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FIG. 15: Scatter plots of the M^^±„ for iV = 11 and 6 = 0.5. Here k = 0.23, 0.245, 0.26, and 0.275 
(running from top-left to bottom-right). 




FIG. 16: Scatter plots of the for the same data set as in Fig. [151 
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in Figs. [13] and [TH At long Monte-Carlo time, when all the M^^ are in the "points" of the 
"stars" , they can be understood semi- quantitatively if 

~ e^™^^"'^°diag(z, i, i, — i, — — i, — — i) x fluctuations , (5.2) 

with n.1 = 77-4 = 1 and n2 = = 2. The order of the diagonal elements in f/^ is unimportant, 
but must be the same for all four links. Note that the matrix in eq. (15.21) does have unit 
determinant as required to be in SU (10), but is traceless and so leads to vanishing P^. 
Fluctuations reduce the magnitude of the M^^, from unity down to about 0.75, and lead 
to the spreading out around the origin. Thus one can understand the behavior for 
these order-parameters as due to the eigenvalues clumping into two subsets (thus breaking 
the symmetry down to (^2)^), and the "locking" of the eigenvalues of the different links 
(breaking the symmetry further down to Z2). Note that the precise form of the eigenvalue 
clumping is dependent on A^. For example, the form in eq. (15. 2p cannot be generalized to 
odd N, for which one eigenvalue is "left out" . This can be used to understand why the 
in the last two panels of Fig. [16] (for which = 11) are not centered on the origin. 

As K is increased further, there is another transition (or transitions) to a phase (or phases) 
in which both the M^j, and the P^ show symmetry breaking (as illustrated by the last panel 
of Figs. [13] and [13]) We have also observed a significant dependence on initial conditions in 
this region. 

In summary, there is a complicated phase structure for k > 0.2, the details of which 
depend on A^. For our purposes, however, the important conclusion is that, in this region, 
the center-symmetry is broken, and so reduction fails. Thus we have not attempted a 
thorough study of this region. 

We close this subsection by comparing our results to those from the 1-loop calculation 
of Ref. [1^. The latter finds that the center-symmetry is broken only for k > 1.4, a much 
larger value than that we find in the single-site model. This large difference may simply be 
due to the difference in the geometries: a single short direction versus four short directions. 



It may also be because Ref. jlj] only estimated energies of simple vacua such as those 
corresponding to an unbroken Zn symmetry or a completely broken Zn symmetry. The 
transition at k ~ 0.04 corresponds to such a breaking, and indeed seems to agree with the 
one-loop result. In contrast to this, the transitions at k > 0.2 involve a more complicated 
breaking of the symmetry, which was not studied in [l^ . 
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FIG. 17: Scatter plots (left panels) and Monte-Carlo time histories (right panels) of 1000 mea- 
surements of Pi (upper panels) and M2,i (lower panels), for N = 10, b = 0.35, and n = 0.155. In 
the right panels, [red] pluses show the real part, while [blue] crosses show the imaginary part. 



F. High statistics study of center-symmetry realization for 0.05 ^ k ^ 0.2. 



In this section we perform more stringent tests to check whether the center-symmetry 
is intact in the physically interesting regime 0.05 ^ k, ^ 0.2. One motivation for doing 
so comes in part from the examples seen in the previous subsection. There we saw that 
simply looking at the Polyakov loops is insufficient because the center-symmetry can be 
only partially broken. Thus it is important to use a set of order parameters sensitive to 
a range of different patterns of symmetry breaking. Another motivation is to push the 
calculation to couplings weaker than b = 0.5, so as to see if there is any barrier to taking 



31 




-1 -0.5 0.5 1 100 200 300 400 500 600 700 800 900 1000 
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FIG. 18: As in Fig. [13 but for b = 1.0 and k = 0.1275. 

the continuum limit in the phase in which reduction holds. In other words, we would like 
to study whether the central "funnel" in Fig. [2] continues up to larger b. Finally, we are 
concerned that the runs discussed so far might be too short to resolve the equilibrium state 
of the theory for some choices of parameters. For example, the time-histories shown in 
Fig. [8] indicate rather long decorrelation times at 6 = 0.5. One would expect this problem 
to worsen as b increases, and indeed we find decorrelation times of 0(50) (as estimated by 
eye) at b = 1.0. We also find qualitative evidence that decorrelation times increase as one 
approaches Kc- 

In order to attempt to address these concerns we performed several long runs, listed 
in Table [TTl consisting of 1000 — 3500 measurements. We begin with a case in which the 
symmetry breaking pattern should be "easy" to resolve based on the results given in previous 
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FIG. 19: As in Fig. [TTl but showing M2,i only, and for b = 1.0, k = 0.09. Upper and lower panels 
show, respectively, results for = 10 (1000 measurements) and N = 13 (3000 measurements — only 
every tenth being shown). 

subsections: the SU{10) theory at 6 = 0.35 and k = 0.1275. This is well within the "funnel" 
and yet not close to Kc. Examples of the scatter plots and time histories are shown in Fig. [TTl 
It appears that these runs are long enough to unambiguously see that Pi and M2,i (and the 
other P^i, and M^,^, for which the plots are similar) are fluctuating around zero. This is 
consistent with the conclusion drawn above, namely that the center-symmetry is intact for 
these values of b and k. 

We now move to weaker coupling. The decorrelation time increases noticeably at 6 = 0.5 
(not shown), although the evidence for the absence of symmetry-breaking remains strong. 
By the time one reaches b = 1.0, however, the results do not have such a clear-cut inter- 
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pretation. This is illustrated in Fig. [181 which shows results at a value of k chosen to be in 
roughly similar relation to as that used in Fig. [T71 so that the quark masses are roughly 
comparable. The scatter plots are not symmetric about the origin, and it is difficult to tell 
from these results alone whether this indicates simply that the run is too short or whether 
the large fluctuations in the time histories are in fact tunneling events between different 
phases in which the symmetry is broken. We think the former possibility more likely, but 
the latter should be kept in mind at this stage. 

One way to differentiate between these two interpretations is to study how the fluctuations 
depend on N. If the symmetry is intact, then (|-P^P) and {\M^^y\^) should vanish as 
as N ^ oo. If, instead, the symmetry is broken, they should tend to a finite value, with 
l/N"^ corrections. We have tried to make this test by comparing SU{10) and SU{13) runs 
at 6 = 1.0, K, = 0.09. An example is shown in Fig. [191 The fluctuations do decrease as 

increases, with (IM^^I) dropping from 0.0850(28) to 0.0620(25). If we rely on these two 
values, then extrapolating in 1/A^^ to = oo yields 0.023(11). This is consistent with no 
symmetry-breaking. 

Another option for studying symmetry-breaking is to study expectation values of the 
operators Kf^, deflned in Sec. IV A[ Since there are a large number (~ 10^) such observables, 
we have to flnd an efficient way to present the results. We proceed by determining the 
signal-to-noise ratio for each n: 

(where AKft is the error in Ka). If the {ZnY symmetry is unbroken, all the expectation 
values should be consistent with zero within errors. Thus we expect to be distributed 
approximately as a Gaussian with width ~ 1. We do not expect an exact Gaussian because 
we are working at flnite A^ and because the observables Kf^ are correlated. Nevertheless, if 
there is symmetry breaking, and some of the observables have non-zero expectation values, 
we expect outliers with |r^| ^ 1. Thus we study many possible realizations of the {Z]\f)^ 
symmetry by looking at the histogram of the r^. Note that since Kfi is a complex number, 
we perform this analysis for both its real and imaginary part, and denote the corresponding 
ratio and histograms by r^ai^mag and -^(rreai.imag)- 

We begin showing in Fig. [20] the histograms if(rreai) for two choices of parameters where 
we know what to expect. These are 
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SU(1 0), b = 0.35, K = 0.1 275, real part SU(1 0). b = 0.5. k = 0.495, real part 




FIG. 20: Histograms of the signal to noise ratio of the real part of Kf^. Left: A case where the 
center symmetry is intact {N = 10, b = 0.35, k = 0.1275). Right: A case where the symmetry is 
broken (A^ = 10, 6 = 0.5, k = 0.495), with the top of the histogram cut off. 

• SU{10) at 6 = 0.35 and k = 0.1275 (1000 measurements), for which all the evidence 
discussed above strongly suggests that the Z^r symmetry is intact (see for example 
Fig.[I7D. 

• SU{10) at b = 0.5 and n = 0.495 (100 measurements), where we have strong evidence 
that the symmetry is broken (see the bottom-right panels of Figs. [13] and [H]). 

For the first choice (left panel) we see the expected Gaussian-like distribution, with almost 
all observables consistent with zero within 2a. By contrast, for the second choice (right 
panel), there are many traces whose signal-to-noise ratio is very large, of O(IOO). This is a 
clear indication that the symmetry is broken (and the pattern of breaking can be deduced 
by determining for which n the signal is significant). 

Now that we have confidence in this method, we apply it to cases where the symmetry 
realization is less clear. Examples of the results are collected in Fig. [2T1 Here the top row 
shows -ff (rj-eai) at 6 = 0.35 and = 10 as we approach closer to Kc (which, from Fig. [TTl 
is at ~ 0.175), and should be compared to the left panel of Fig. [201 The middle row of 
Fig. [21] shows if(ri.eai) at weaker coupling (6 = 1.0, still A^ = 10) both away from (k = 0.09) 
and close to {k, = 0.1275) Hc- These are from the same data sets as those illustrated in 
the upper panels of Fig. [191 and all panels of Fig. [18], respectively. Finally, the bottom row 
shows if(rreai) and -ff (nmag) (left and right panels) at 6 = 1.0, k, = 0.09, but now at A^ = 13 
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SU(10), b = 0.35, K = 0.15, real part 



SU(1 0), b = 0.35, K = 0.1 55, real part 




-1 1 



SU(10),b = 1,K = 0.09, real part 




SU(13),b = 1,K = 0.09, real part 
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600 - 



SU(10),b = 1,K = 0.1275, real part 




SU(1 3), b = 1 , K = 0.09, imaginary part 




FIG. 21: Histograms of for various parameters of interest. From top-left to bottom-right, 
parameters are (i) = 10, 6 = 0.35, n = 0.150 (real part), (ii) = 10, 6 = 0.35, n = 0.155 (real 
part), (iii) N = m,h = 1.0, k = 0.09 (real part), (iv) iV = 10, 6 = 1.0, k = 0.1275 (real part), (v) 
AT = 13^ 6 = 1.0, K = 0.09 (real part), (vi) iV = 13, 6 = 1.0, k = 0.09 (imaginary part). 
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(corresponding to the data in the lower panels of Fig. [T9|) . 

In none of these cases is there any evidence for outliers, and we conclude that it is unlikely 
that the center-symmetry breaks in "funnel" region, at least up to 6 = 1.0. We note that 
this histogram method appears to be a more powerful tool than looking at individual scatter 
plots and time-histories. 

VI. RESULTS OF PHYSICAL INTEREST 

Having found evidence that reduction holds in the interesting region on either side of the 
putative critical k, we now make a first attempt at extracting quantities of physical interest. 
As — s> oo, the results we find should hold also for the large- infinite- volume gauge theory 
with one adjoint quark (as long as one uses the same action and the same values of b and 
k). As k, approaches k^, the quarks become light, and their contribution to the dynamics 
becomes important. Conversely, as k moves away from Kc (towards zero, say), the quarks 
become heavy (compared to its dynamically generated scale) and the dynamics approaches 
that of the pure-gauge theory. The fact that reduction appears to hold down to small values, 
K ~ 0.05, where we expect the quarks to be very heavy, indicates that we may well be able 
to use reduction to study the pure-gauge theory using this "adjoint deformation" . 

What we would like to do is use known results from the infinite-volume pure-gauge 
theory to see how close we are to N = oo.^^ Unfortunately, we cannot do any quantitative 
comparisons at this stage. For one thing, since we are at this stage unable to measure pion 
masses, we do not know where the boundary between light and heavy quarks lies. And even 
if we had determined that the quarks are heavy, so that the long-distance physics was that 
of the pure-gauge theory, the presence of a non-zero k would lead to additional terms in the 
gauge action, including, for example, the trace of the plaquette in the adjoint representation. 
In other words, reduction would match our single-site model to a pure-gauge theory with a 
different gauge action, with the additional terms entering at 0{k^). 

Despite these drawbacks, we think it useful to attempt a large- A^ extrapolation for some 
quantities in order to make a qualitative comparison to the pure-gauge theory. We do so for 
the average plaquette and for the distribution of eigenvalues of a quenched overlap fermion 

We are not aware of any results for the large-volume gauge theory with Nf — 1 dynamical adjoint quarks 
with which to compare. 
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in the fundamental representation. We also attempted to extract the string tension from 
the e"^*^ dependence of Wilson-loops, with A the area. (The loops are calculated using 
the reduction prescription of Ref. 2|). We find that the Wilson-loop expectation values do 
drop as A increases, but only for a short window, A < A^, after which they start to grow. 
This growth is presumably due to corrections. Although we find that the upper edge 
of the window, Ac, grows with N, the window is too small at our values of to extract 
a string-tension. This problem could be resolved either by using larger values of or by 
developing variational techniques to extract the tension at short distances. We leave this 
for future studies. 



A. Average plaquette 

We begin by comparing the values of the plaquette. We focus on three (6, k) values 
at which we have good statistical control for A^ = 8 — 15: (0.5,0.09), (0.5,0.1275) and 
(1,0.09). The results are collected in Table UTTl We also include in the Table the value for 



the pure-gauge theory in the large A^ limit. For 6 = 0.5 this is given by 



44| 
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«(k = 0,6 = 0.5) -0.7182, (6.1) 

while for larger values of b, we use three-loop perturbation theory (taken from, for example, 
Ref. l22j) 

, ^ b^oo , 1 0.653687 0.4066406 , ,^ 

MK = 0,O — > 1 — — h... 6.2 

^ ' ^ 8 6 128 62 512 63 v ; 

= 0.8692... at 6 = 1.0. (6.3) 

At 6 = 1.0 the 1/6'^ term contributes about 0.1% to the plaquette value and so we estimate 
the error in the result to be ~ 0.0008. 

What we see is that the results at k = 0.09 (for both values of 6) approach the pure-gauge 
theory value from above as A^ — * oo, while that at k = 0.1275 (where the quark mass is 
lighter) approaches from below. To see how close they approach this value, and to study 
the nature of the 1/A^ expansion, we fit our data to 

u{N) = u{oo) + — , with g = 1 or 2. (6.4) 
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The error on this number is very smaU and we can safely assume that it is zero in the discussion below. 
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b 


u 




SU{8) 


SU{W) 


SU{11) 


SU{13) 


SU{15) 


0.5 


0.718 


0.09 


0.7460(20) 


0.7429(20) 


0.7420(12) 


0.7388(22) 


0.7362(12) 


0.5 


0.718 


0.1275 


0.6836(12) 


0.6877(7) 


0.6959(18) 


0.6974(18) 


0.6992(16) 


1 


0.870 


0.09 


0.88304(10) 


0.8700(4) 


0.8798(14) 


0.8779(6) 


0.8774(8) 



TABLE III: Comparison of plaquette expectation values u of large- pure-gauge theory with those 
obtained in our single-site simulations. 



The 1/A^^ fit is appropriate if we are in the asymptotic regime, while the 1/A^ fit is an attempt 
to mock up the behavior if we are far from asymptotia. The resulting extrapolations are 
shown in Fig. [22], with fit parameters given in Table IIVI We observe that the linear and 
quadratic fits are of comparable quality, indicating that we need a wider range of to pin 
down the appropriate fitting form. One can perhaps use the difference between these fits as 
a crude estimate of the extrapolation uncertainty. The rather large values in the second 
row of the table reflect the scatter of our data around the fit lines (see Fig. I22l) and might 
indicate underestimated statistical errors or a competition between multiple terms in the 
1/A^^ expansion. 

As noted above, reduction does not predict that the extrapolated single-site plaquette 
values should agree with those in the pure-gauge theory at infinite-volume, only that they 
should be close for small k. It is therefore slightly surprising that we find better agreement 
for K = 0.1275 than for 0.09 at 6 = 0.5. Further work will be required to determine the 
significance of this finding. 



Data set 




Type of fit 


u{oo) 


A 


xVd.o.f. 


b = 0.5, K = 0.09, u = 


0.718 


Linear 

Quadratic 


0.7255(32) 
0.7332(17) 


0.171(37) 
0.91(20) 


1.45/1 
1.1/1 


b = 0.5, K = 0.1275, u -- 


= 0.718 


Linear 
Quadratic 


0.7183(30) 
0.7085(16) 


-0.291(30) 
-1.43(15) 


12/3 
16/3 


b = 0.5, K = 0.09, u = 


0.870 


Linear 
Quadratic 


0.8700(12) 
0.8744(5) 


0.101(16) 
0.559(77) 


0.97/3 
0.64/3 



TABLE IV: Results of extrapolations of plaquette to A^ = oo, obtained from fitting the data in 
Table Im] to the form Eq. (lOl) . 
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kappa = 0.09, b = 0.5 
0.755 1 — , , , 



0.75- 




0.67' — ' • < • 1 

0.05 0.1 0.15 0.2 

1/N 

kappa = 0.09, b = 1 
0.888 1 — • • • • 1 



0.886 - 




FIG. 22: Large-A^ extrapolations of the average plaquette for (6, k) = (0.5,0.09), (0.5,0.1275) and 
(1.0,0.09) (running from top to bottom). The filled circle [red] at = oo is u, the value in the 
large- A'^ limit of the infinite- volume pure-gauge theory. 
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B. Dirac spectrum of fundamental fermions 



One drawback of the average plaquette is that its value is dominated by short-distance 
physics (i.e. gauge fluctuations with wavelenths ~ We consider in this section a 

quantity that is sensitive to long-distance physics, and thus serves as a better test of whether 
our values of are large enough to extract long-distance quantities. 



The idea, proposed in Ref. 



45| , is to probe the large- theory using the eigenvalue spec 



trum of valence fermions in the fundamental representation. As discussed in Ref. [45|, it 
is legitimate in this context to quench fundamental representation fermions in the large- A^ 
limit. The specific proposal is to calculate the distribution of the low-lying eigenvalues 
of the quenched overlap Dirac operator and compare them to the predictions of random 
matrix theory (RMT). This constitutes a test that the chiral-symmetry-breaking dynamics 
of large- A^ QCD are being correctly reproduced, because the RMT predictions can be de- 
rived from QCD if the eigenvalues are in the so-called "epsilon-regime" and if some other 
conditions hold [471]. One can furthermore extract a value for the condensate, {qq)/N, from 
this comparison. 



This approach is used in Ref. [45| in the context of partial volume reduction, in which 
one simulates SU{N) pure-gauge theories in boxes of physical size of 0(1 fm), which are 
found to be large enough to satisfy volume-independence [l8|. It is argued in Ref. j45| that 
the eigenvalue densities of the valence Dirac operator of the (partially) reduced theory are 
legitimate quantities to be compared with those of RMT as long as A^ is large enough. 
Specifically, one should expect that the smallest eigenvalues are described by RMT once 
NL'^ is large enough. This can be achieved either by increasing L or by increasing A^. 



A particularly useful quantity considered in Ref. [45|] was the distribution, P{p), of the 
ratio between the first and second eigenvalues of the overlap Dirac operator. This has the 
advantage, compared to the distributions of individual eigenvalues, that the RMT prediction 
is parameter-free (i.e. is independent of the value of the condensate). Thus it can be used 
as a gauge of whether NL'^ is large enough. Indeed, using this quantity, Ref. j45| found that 
the measured -P(p) agrees with the RMT prediction on an L = 6 lattice only for A^ > 23. 
In our case we have L = 1 and so we probably need even larger values of A^ in order to 
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See Ref. [46j for further discussion of the conditions under which such quenching is, and is not, justified. 
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FIG. 23: The quantity P{p) (see text) for SU{1^) and b = 0.35 compared to RMT (solid curves). 
The values of k are 0.1275 (left panel) and 0.155 (right panel). 

see -P(p) asymptote to the its RMT form. Note that this expectation is justified only for 
values of n that correspond to adjoint fermions heavier than the dynamical scale of the gauge 
theory. If the adjoint fermions become light, then their determinant will alter the expected 
distribution. 



We thus calculate P{p) for valence overlap fermions (using the conventions of Ref. |45| . 
and taking Mq = —1.5). Since the dimension of the Dirac matrix of the fundamental fermions 
is modest, 36 — 60, we can construct it exactly, without approximating the sign function 
involved in its definition. The parameters for which we calculated P{p) were N = 10, 
b = 0.35 and k = 0.1275,0.155. These we expect to correspond to relatively heavy and 
moderately light fermions, respectively, based solely on their proximity to ~ 0.175. We 
used 1000 gauge configurations, all of which we found to have zero topological charge (using 
the index theorem). We present our results for P{p) in Fi g. [23l together with a solid curve 
that is the analytic formula of RMT reproduced from Ref. 451]. 



The disagreement between our data and the RMT prediction is clear. The form of the 



disagreement is, in fact, similar to that seen in Ref. [45|| for small values of A^. This leads 
us to conclude our values of are too small for the lowest two eigenvalues to be in the 
epsilon-regime. It would be of considerable interest to extend this calculation to larger A^. 
Another obvious step is to perform a comparable study with a Dirac operator of a fermion 
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We found that a few configurations had IQtopl = 1 if we changed Mq. 
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in the adjoint representation. A straight-forward generahzation of the arguments in Ref. 45 1 
shows that to be in the epsilon-regime now requires N'^L'^ to be large enough, which is easier 
to satisfy, and it may be that our modest values of suffice. One complication is that our 
fermion action, which uses the Wilson-Dirac operator, does not preserve chiral symmetry, so 
that we cannot making a direct comparison with RMT unless we include lattice artefacts. 



It may be possible to use the approach of Ref. [48|, however, to study the condensate. It 



also may be possible to use a valence overlap operator for the adjoint fermions. 

VII. SUMMARY, CONCLUSIONS AND FUTURE DIRECTIONS 

In this paper we have taken a step towards exploring large- A^ QCD with fermions in two- 
index representations. In these theories the number of both gauge and fermionic degrees of 
freedom grows as 0{N'^), so that the latter contribute to the dynamics even when A^ — > oo. 
We do not study these theories directly, but rather use large- A^ equivalences to relate them 
to a much simpler theory, QCD with adjoint fermions defined on a single site. These equiv- 



alences follow from a combination of orbifold and orientifold projections [g]. Specifically, we 
choose to work with a single Dirac adjoint fermion which, through these projections, corre- 
sponds to a gauge theory with two Dirac fermions in the antisymmetric representation. For 
A^ = 3 the latter theory becomes physical (3-color) QCD with two degenerate Dirac fermions 
in the fundamental representation, and this makes our study relevant phenomenologically. 

Clearly, for our approach to work it is necessary that the above-noted equivalences hold. 
This in turn requires that the ground state of the reduced theory is symmetric under the 
{ZnY center-symmetry of the theory. The present paper is focused on determining, using 
Monte-Carlo simulations, the regime in the parameter space of the single-site theory within 
which the center symmetry remains unbroken. Our simulations were performed with 8 < 
A^ < 15 at a variety of lattice spacings and quark masses. The observables we measure 
include order parameters for the breaking of {ZjsiY symmetry, and in some instances we 
gather large data samples, allowing the calculation of ~ 10^ different order parameters that 
probe many potential patterns of symmetry-breaking. 

We find strong evidence that the center symmetry is intact in an extended region of the 
lattice parameter space, a region that includes the "critical" line along which we expect 
that the fermions have their minimum mass. Our results for the phase-diagram depend very 
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weakly on N, suggesting that they apply also when N oo. In particular, the relatively 
small values of that we use appear large enough to observe the first-order transition line at 
Kc, despite the fact that this becomes a true transition only when oo. Our results are 
consistent with the region of unbroken center-symmetry extending towards the continuum 
limit, so that the phase diagram is consistent with that conjectured in Fig. [2], although we 
cannot rule out that this region shrinks when 6 > 1. 

An important finding is that the center symmetry does not break until the physical 
fermion mass becomes very heavy, likely at the cut-off scale (in approximate agreement with 
the analytic estimate Ref. jl^). For example, at 6 = 0.5, where Kc ~ 0.15, the symmetry is 
unbroken for k = 0.06 (see Figs. E] and [7]), so m = (Zm/a)[l/(2/t) — 1/(2kc)] ~ (Z^/a) x 5. 
Since we expect ~ 0{1), the fermion mass is of 0(1/ a). Thus there appears to be an 
overlap of the region in which reduction holds and that in which the long-distance physics 
of the corresponding large- volume theory is that of large- A^ pure-gauge theory. This finding 
opens a window to the study of the pure-gauge theory using reduction, which was the original 
idea behind the proposal of Eguchi and Kawai [2]. It seems likely that what is happening 
here is that the heavy fermions would, if integrated out, induce a tower of interactions 
between Polyakov loops that is similar to the tower of double-trace interactions proposed in 
Ref. [12i] to stabilize the center symmetry. 

As already noted, our evidence for the absence of center-symmetry breaking becomes less 
strong at the smallest coupling we consider, b = 1.0. This is an extremely small coupling, 
corresponding to /3 = 18 if A^ = 3. It is much smaller than the values for which we envision 
performing useful measurements of physical observables {b ^ 0.35). The issues that arise 
at 6 = 1.0 are that there are large fluctuations in Monte-Carlo time histories, and long 
auto-correlation times, making it hard to unambiguously determine the equilibrium state. 
It is for these couplings that the use of the ~ 10^ order parameters becomes particularly 
useful, allowing us to try and tease out evidence of symmetry-breaking. We find none. 

A general lesson we have learned is that, when shrinking more than one Euclidean di- 
rection, the center-symmetry sometimes breaks in quite nontrivial ways. For example, in 
some parts of the lattice phase diagram we observed ground states for which the expectation 
values of the Polyakov loop are consistent with zero, while other order parameters, which 
measure correlations between different euclidean directions, have non-zero averages. This is 
similar to the behavior we observed in the quenched Eguchi-Kawai model [2j], and shows 
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that it is insufficient to measure only Polyakov loops (or any power thereof) when studying 
center-symmetry breaking. 

Of course, finding that reduction holds is only the ffist step. Our ultimate aim is to use 
the single-site theory as a tool for learning about physical quantities of large- QCD with 
one flavor of adjoint fermions, and large- QCD with two flavors of fermions in the antisym- 
metric representation. Reduction allows one to calculate expectation values of Wilson loops 
and connected correlation functions of certain fermionic operators, and thus to determine 
the string tension and certain "hadron" masses, as well as glueball-gg mixing, etc.. From 
a practical point of view, however, the key question is this: What value of is needed to 
obtain results with controlled corrections? We have made a first step at answering this 
question by looking at two variables — the average plaquette and the eigenvalue densities of 
a fundamental-representation massless fermion. The former is not itself a physical quantity, 
as it is dominated by ultraviolet fluctuations, but it is simple to calculate and can give an 
indication of the behavior. While our calculations using = 8 — 15 cannot pin down 
the form of the 1/A^ dependence, it does appear that the results at A^ = 15 lie within a few 
percent of those at A^ = oo. 

The eigenvalue densities provide a test of whether the single-site gauge configurations 
can reproduce the long-distance physics of chiral-symmetry breaking. In particular, in their 
"epsilon- regime" , the low lying eigenvalues of the Dirac operator have correlated distribu- 
tions that are predicted by random-matrix theory (RMT). We find distributions which differ 



significantly from the RMT predictions, and, based on the results of Ref. j45|, take this as 
evidence that our values of A^ are too small to probe the epsilon-regime. 

It is important to stress that, although we find that A^ = 8 — 15 are too small for this 
particular observable, these values of A^ do appear to be sufficient to study the nature of 
the phase diagram, which is the main goal of this paper. 

Looking forward, the prospects for using reduction as a quantitative numerical tool clearly 
depend on developing or implementing algorithms with a less formidable A^ dependence than 
the A^^ scaling of our method. This will allow us to simulate much larger values of A^, and 
see if the phase diagram presented in this paper for A^ < 15 is indeed indicative of its A^ = cxo 
limit. We are considering implementing a hybrid-Monte- Carlo. It also is clear that to move 
forward one needs greater computing resources. In that regard, one issue to be faced is the 
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lack of any obvious way of parallelizing the code. Here the advantage of reduction — packing 
as much information as possible into the link matrices — leads to a computational problem. 

A less ambitious direction of further study is to make more extensive and systematic 
measurements of observables on the lattices we already have at hand. We have in mind in 
particular studying the correlators of hadrons composed of adjoint representation fermions. 
This would allow us to check our hypothesis that the "critical" line corresponds to the 
minimal mass of the pions. We would also like to measure the eigenvalue densities for 
adjoint-representation fermions. These are expected to enter their epsilon-regime for values 
of that are parametrically smaller than those required in the fundamental fermions case. 

Finally, we recall that using reduction is also of interest for other values of the number 
of Dirac flavors. In particular, for Nf = 1/2, the chiral limit of the theory is related to the 
A/" = 1 SUSY gauge-theory, while, for Nf = 2, the theory is perhaps conformal. If reduction 
holds for these cases as well, then it provides a method to study the large- A^ limits of these 
theories. 



Acknowledgments 

We thank H. Neuberger for bringing Ref. Sj to our attention, and Joyce Myers for dis- 
cussions. This study was supported in part by the U.S. Department of Energy under Grant 
No. DE-FG02-96ER40956. 



[1] M. Teper, PoS LATTICE2008, 022 (2008) arXiv:0812.0"n85 ' [hep-lat]. 

[2] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 1063 (1982). 

[3] H. Neuberger, Annales Henri Poincare 4, S147 (2003) |arXiv:hep-th/021 2097 1 . 

[4] P. Kovtun, M. Unsal and L. G. Yaffe, JHEP 0507, 008 (2005) | arXiv:hep-th/0411177j . 

[5] A. Armoni, M. Shifman and G. Veneziano, Phys. Rev. Lett. 91, 191601 (2003) 

|arXiv:hep-th/0307097] ; Nucl. Phys. B 667, 170 (2003) | arXiv:hep-th/030216 3|; Phys. Rev. D 

71, 045015 (2005) |arXiv:hep-th/0412203] . 
[6] P. Kovtun, M. Unsal and L. G. Yaffe, JHEP 0706, 019 (2007) |a rXiv:hep-th/070 202l]. 
[7] M. Unsal and L. G. Yaffe, Phys. Rev. D 74, 105019 (2006) | arXiv:hep-th/0608180] . 



46 



[8] G. Bhanot, U. M. Heller and H. Neuberger, Phys. Lett. B 113, 47 (1982). 
[9] V. A. Kazakov and A. A. Migdal, Phys. Lett. B 116, 423 (1982). 
[10] M. Okawa, Phys. Rev. Lett. 49, 353 (1982). 

[11] P. N. Meisinger and M. C. Ogilvie, Phys. Rev. D 65, 056013 (2002) | arXiv:hep-ph/0108026] ; 

J. C. Myers and M. C. Ogilvie, arXiv: 0903. 46381 [heo-th]: 

T. J. Hollowood and J. C. Mvers. [aFXw:0907.3665l [heo-th]. 
[12] M. Unsal and L. G. Yaffe, Phys. Rev. D 78, 065035 (2008) [a rXiv: 0803. 0344] [hep-th]]. 
[13] J. C. Myers and M. C. Ogilvie, Phys. Rev. D 77, 125030 (2008) tarXiv:0707.1869 [hep-lat]]. 
[14] B. Bringoltz. larXiv:0905.240 6' [hep-lat]. 

[15] P. F. Bedaque, M. L BuchofF, A. Cherman and R. P. Springer. [arXiv:0904.0277] [hep-th]. 
[16] J. Ambjorn, K. N. Anagnostopoulos, W. Bietenholz, T. Hotta and J. Nishimura, JHEP 0007, 

013 (2000) |arX iv:hep-th /0003208]. 
[17] M. Hanada, L. Mannelh and Y. Matsuo. larXivi0905.2995 [hep-th]. 

[18] J. Kiskis, R. Narayanan and H. Neuberger, Phys. Rev. D 66, 025019 (2002) 
|arXiv:hep-la t /0203()0F| . 

[19] A. Gonzalez-Arroyo, R. Narayanan and H. Neuberger, Phys. Lett. B 631, 133 (2005) 

|arXiv:hep-lat/05 09074] . 
[20] A. Gonzalez-Arroyo and M. Okawa, Phys. Lett. B 120, 174 (1983). 

[21] W. Bietenholz, J. Nishimura, Y. Susaki and J. Volkholz, JHEP 0610, 042 (2006) 

|arXiv:hep-t h70608072 1 . 
[22] M. Teper and H. Vairinhos, Phys. Lett. B 652, 359 (2007) [arXiv:hep-th/0612097| . 
[23] T. Azeyanagi, M. Hanada, T. Hirata and T. Ishikawa, JHEP 0801, 025 (2008) 

|arXiv:071 1.1925 [hep-lat]]. 
[24] B. Bringoltz and S. R. Sharpe, Phys. Rev. D 78, 034507 (2008) [ arXiv:0805.2T46] [hep-lat]]. 
[25] J. Giedt, R. Brower, S. Catterall, G. T. Fleming and P. Vranas, Phys. Rev. D 79, 025015 

(2009) [arXiv:0810.5746l [hep-lat]]. 
[26] M. G. Endres. [arXiv:0902.4267 ] [hep-lat]. 
[27] J. Giedt. [aIXw:0903.2443l [hep-lat]. 

[28] L. Del Debbio, A. Patella and C. Pica. 'a rXiv:0805.2058l [hep-lat]: 

A. J. Hietanen, K. Rummukainen and K. Tuominen. larXiv:0904. 08641 [hep-lat]: 
G. Cossu and M. D'Eha, larXiv:0904. 13531 [hep-lat]. 



47 



[29] S. Catterall, J. Giedt, F. Sannino and J. Schneible, JHEP 0811, 009 (2008) iarXiv:0807.0792] 
[hep-lat]]. 

[30] M. E. Peskin, Nucl. Phys. B 175, 197 (1980). 

[31] H. Leutwyler and A. V. Smilga, Phys. Rev. D 46, 5607 (1992). 

[32] J. B. Kogut, J. Polonyi, H. W. Wyld and D. K. Sinclair, Phys. Rev. Lett. 54, 1980 (1985). 
[33] F. Karsch and M. Lutgemeier, Nucl. Phys. B 550, 449 (1999) |arXiA^p-lat /9812 023 1 
[34] I. Montvay, Int. J. Mod. Phys. A 17, 2377 (2002) [arXiv:he p-lat7011 2007 1 . 
[35] K. Demmouche, F. Farchioni, A. Ferling, G. Munster, J. Wuilloud, I. Montvay and 

E. E. Scholz, larXiv:0810.0T44] [hep-lat]. 
[36] S. R. Sharpe and R. L. Singleton, Jr., Phys. Rev. D 58, 074501 (1998) jarXiv:hep-lat/98 04028 j . 
[37] L. Del Debbio, M. T. Frandsen, H. Panagopoulos and F. Sannino, JHEP 0806, 007 (2008) 

|arXiv:0802.0"89n [hep-lat]] . 
[38] S. Aoki, Phys. Rev. D 30, 2653 (1984). 
[39] S. Aoki, Phys. Rev. D 33, 2399 (1986). 
[40] M. Teper, private communication. 

[41] M. Campostrini, Nucl. Phys. Proc. Suppl. 73, 724 (1999) [arXiv:hep-lat/9809072| . 

[42] J. Kiskis, R. Narayanan and H. Neuberger, Phys. Lett. B 574, 65 (2003) 

|arXiv:hep-Iat/0308033] . 
[43] N. Cabibbo and E. Marinari, Phys. Lett. B 119 (1982) 387. 
[44] M. Teper and H. Vairinhos, private communications, (2008). 

[45] R. Narayanan and H. Neuberger, Nucl. Phys. B 696, 107 (2004) |arXiv:hep-lat /0405025] . 
[46] B. Bringoltz and L. G. YafFe, in progress. 

[47] J. C. Osborn, D. Toublan and J. J. M. Verbaarschot, Nucl. Phys. B 540, 317 (1999) 
|arXiv:hep-th/9806lT0] ; 

P. H. Damgaard, J. C. Osborn, D. Toublan and J. J. M. Verbaarschot, Nucl. Phys. B 547, 
305 (1999) | arXiv:hep-th79811212] . 
[48] L. Giusti and M. Luscher, JHEP 0903, 013 (2009) |arXiv:0812.3638] [hep-lat]]. 



48 



